New analytic unitarization schemes 
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We consider two well-known classes of unitarization of Born amplitudes of hadron elastic scat- 
tering. The standard class, which saturates at the black disk limit includes the standard eikonal 
representation, while the other class, which goes beyond the black disk limit to reach the full 
unitarity circle, includes the U matrix. It is shown that the basic properties of these schemes are in- 
dependent of the functional form used for the unitarisation, and that U-matrix and eikonal schemes 
can be extended to have similar properties. A common form of unitarization is proposed interpo- 
lating between both classes. The correspondence with different nonlinear equations are also briefly 
examined. 



I. INTRODUCTION 

At low energies, hadron scattering can be described 
by one-Reggeon exchange terms. But as the Pomeron 
term(s) grow(s) with energy, these exchanges will even- 
tually violate unitarity. To see this, one can switch to 
partial waves [![, or to the representation in impact pa- 
rameter b. As s grows, one needs to sum many partial 
waves with I ~ |b|^/i — > oo. The summation over I then 
becomes an integration over b. 

The partial wave G(s, b) has two regimes. First of all, 
it can reach maximum inelasticity. In this case, G(s, b) = 
1, and half of the interactions are inelastic. The center of 
the protons then becomes black, and multiple exchanges, 
i.e. cuts in the complex J plane, become important. 
Second, the partial wave can later reach the full unitarity 
limit G(s, b) = 2. 

The maximum inelasticity limit may be reached in 
or pp scattering a little above the Tevatron energy 
so that one expects cuts to be important in the descrip- 
tion of soft interactions at the CERN LHC. The inclu- 
sion of these goes under the name of unitarization. It is a 
formidable task to calculate the contribution of cuts, as 
not only multiple Pomeron exchanges must be calculated, 
but also multiple Pomeron vertices. 

Different schemes have been proposed, and we want in 
this paper to show that the general properties of the am- 
plitude do not heavily depend on the scheme, but rather 
on what assumes for the inelastic contribution at high 
energy. We shall limit ourselves to two popular schemes: 
the eikonal and the U matrix, and show that simple ex- 
tensions of each lead to similar properties. 

In Sec. I, we remind the reader of the simple require- 
ments coming from unitarity, and examine in Sec. II the 
two schemes. In Sec. Ill, we show that it is possible to 
obtain the properties of the eikonal by extending the U - 
matrix scheme, whereas in Sec. IV, we show the reverse, 
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i.e. that one can extend the eikonal to mimic the behav- 
ior of the U matrix, at least for physical amplitudes. 



II. UNITARITY 

At high energy, we can start with the elastic scatter- 
ing amplitude a(s,t), related to the elastic cross section 
though 



da 



1 



dt Wirs 2 1 
One can then Fourier transform a to b space 

[ ■ > J (2tt)2 2s 
which leads to the expressions 



citot = 2 J cPbImG(s,b), 

(Tel = J d 2 b |G( S ,b)| 2 , 



(1) 

(2) 

(3) 
(4) 



where we have assumed that the spin-flip contribution to 
the elastic cross section is negligible. One can then write 
the square of the S'-matrix density S(s, b) = 1 + iG(s, b) 
as 



\S(s, b)| 2 = 1 - 21m G(s, b) + \G(s, b)|' 



(5) 



Unitarity demands that | aS' ( t») | 2 1 < 1, the difference com- 
ing from inelastic channels: 



7? m ( S ,b) = l-|S(s,b)| 2 >0. 



(0) 



There are several ways to represent the unit circle. 
First of all, one can map the upper complex plane into a 
circle via a complex exponential 



S(s,h) = exp(iz(s,b)) with Imz(s,b) > 0, 



(7) 



This maps in fact an infinite number of strips with 2n7r < 
Rez(s, b) < 2(n + 1)7T each onto the unit circle. 
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It is also possible to use a one-to-one map through a 
Mobius transform and write 

S(a, b) = \ + _\ Z ^ h ^y with W(s, b) > 0. (8) 

Other representations are possible, but we shall concen- 
trate on these two in the following section. 



III. UNITARIZATION 

The physical amplitude lies within the unitarity circle, 
so that the associated S matrix can always be represented 
by Eqs. {7} and ©. The unitarization scheme comes in 
once one identifies z or z' with the one-Reggeon exchange 
amplitude. One then considers ([7]) and (J8]) as series ex- 
pansions in n-Reggeon exchanges, so that their first term 
must give 1 + «x( s j 

Indeed, if one writes the one-Reggeon exchange ampli- 
tude as x(s,b), then assuming z = x m ® leads to the 
well-known eikonal representation: 



G(s, b) = i(l - exp(ix(s, b)). 



(9) 



This scheme can be derived in QED and other field the- 
ories 0, 0, U @ or m potential theory 0]. It can be 
extended to include diffractive channels Q. It leads at 
asymptotic energies (s — > oo) to the limit cr e i/<T- me i = 1, 
i.e. to maximum inelasticity. 

The other unitarization scheme considered here is the 
[/-matrix representation 0, liol . [llj where one identifies z' 
in with x(s, b)/2, to match the one-Reggeon exchange 



G(s,b) = 



X(s,b) 



ix{s,b)/2 



(10) 



In this scheme, S^Sjb) tends to —1 when s — > oo and b 
is finite, so that the inelastic partial wave r?i n (s,b) tends 
to 0: the ratio a e \/a m c\ vanishes asymptotically. 

Both schemes have the same development at second 
order in %, and differ only in the rest of the series. 

It must be noted however that the resummation must 
lead to an amplitude within the unitarity circle, but there 
is no reason to assume that it maps the entire complex 
plane to the circle. Hence, one can easily extend both 
schemes through a change in the strength of successive 
scattering. This gives the extended eikonal schemes [H, 
li 



G(s,b) = -(1 - exp(iujx(s,b)) 

and the extended [/-matrix schemes 

X(s,b) 



G(s,b) = 



1 -iu)'x(s,b)' 



(11) 



(12) 



It is straightforward to check that using oj > 1 or oj' > 
1/2 maps any amplitude x into the unitarity circle. 

We shall now show that the various possibilities can be 
grouped into two wide classes of unitarization schemes, 
and that the exact form matters little. 



IV. SHADOWING 

As we have seen, the eikonal predicts that at high en- 
ergy the inelastic component of the cross section will be 
maximal, T]- m = 1. This in turn leads to |5(oo,b)| = 
and G(oo, b) = i. To reach this regime via an extended 
U matrix, one needs to choose u>' = 1 in (fT2|) . 

The inelasticity will then be 



mn(s,b) 



2Im X ( g ,b) + | X (s,b)| 2 
l + 2Imx(s,b) + | X (s,b) 



(13) 



1 



It can easily be seen that s — > oo leads to <7i ne i/cei 
and crci/ctot — * 1/2. 

Differently stated, this extended U matrix representa- 
tion has the standard black disk limit. 
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FIG. 1: Inelasticity for the eikonal and extended U matrix at 
Y^s = 1.8 TeV (upper figure) and sfs = 14 TeV (lower figure). 

Throughout this paper, we shall use as an example of 
one-Reggeon exchange amplitude a hard Pomeron term, 
with a parametrization 



1500 GeV 



0.9 



exp{-b 2 /(9 GeV' 2 ) 



x I (i + tan 



0.45tt 



(14) 



This amplitude reaches the black disk limit at 1500 GeV, 
similarly to the model of ref . [l5| , and has a dependence 
in t similar to that of pp scattering. We also neglect the 
effect of shrinkage, which is small for a hard Pomeron. 
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In Fig. 1, one can compare the inelasticity in the case 
of the eikonal and in that of the extended U matrix. One 
clearly sees that the generic features of both schemes are 
very close. 

The extended U matrix and the eikonal scheme are 
different representations of a wider class of unitarization 
procedures with a standard black disk limit. Indeed, we 
can extend (TT21 to 



G(*,b)=i[l 



1 



(1 - i x ( Sl b)/l) 7 



(15) 



If 7 = 1 this form leads to the extended U matrix while, 
for 7 — » ao, we obtain the standard eikonal. When 7 
varies from 1 to 00, we obtain different forms of unita- 
rization which all lead to a black disk limit, and the am- 
plitude G does not change anymore once it has reached 
its maximum value. In Fig. 2, we show the inelasticity 
r](s, b) for different values of 7 , again in the case of the 
hard Pomeron input of Eq. (CO 
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FIG. 2: Inelasticity rji n (s, b) for various values of 7, for y/s : 
1.8 TeV (upper line) and y/a = 14 TeV (lower line). 



can also 
* 0, as in 



It may be worth pointing out that Eq. ((To 
lead to a noninteracting theory in the limit n ■ 
this case G(s, b) — > 0. 



V. FULL UNITARITY LIMIT 



The standard U -matrix scheme (ITOl) was intensively ex- 
plored in [Tlj in the partial- wave language. In the impact 
parameter representation, the properties of the U -matrix 
are explored in . 



For a purely imaginary one-Reggeon exchange, when 
Imx(s, b) goes from to 00, the S'-matrix varies in the 
interval [—1, +1], and the amplitude G goes from to 2i. 
Hence in this case the full unitarity limit can be reached. 

This form of unitarization leads to unusual properties 
at super- high energies as was shown in [l6l |. In this rep- 
resentation 



Cei(s) 







2 




2 





1 -i 



and 



Cinel(s) 



Im 



2 



X(s,b) 
2 



db 



(16) 



1 



X(*,b) 
2 



db. (17) 







so that, when Imx(s,b) — > 00, one gets (Tinci/ctot 
and er c i/f7tot — > 1- 

It is often considered [TjJ that these properties are in- 
trinsic to the Mobius projection of \ onto the unitarity 
circle. We want to show now that, in fact, extended 
eikonals can lead to the same properties for the unita- 
rized amplitude. 

We have seen that choosing ui > 1 in (fTTj) guarantees 
that any amplitude would be unitarized. However, one 
must be concerned with the physical amplitude, and it 
is not needed to map the whole complex plane into the 
unitarity circle. This means that for some specific choices 
of one-reggeon exchange x, one can extend the range of 
values of to, and restrict oneself to part of the complex \ 
plane. Unitarity, in this case, leads to the condition 



cos(ojRex) > 



(2w- l)e 



Lolmx 



2(1 -w) 



(18) 



So we see that for u < 1/2, the second term of the nu- 
merator will guarantee the inequality for sufficiently large 
Imx- 

We show in Fig. [3] the region allowed in the case 
lu = 0.525 together with a curve showing the amplitude 
corresponding to the exchange of one hard pomeron with 
intercept 1.45. This is of course an extreme curve, corre- 
sponding to a ratio Rex/Imx of 0.73. Any physical am- 
plitude will include softer intercepts, and will lie above 
the hard pomeron line. So we see that, in practice, eikon- 
als can be extended to values of oj between 1/2 and 1. 

But at high energy, such eikonals have all the basic 
properties of the [/-matrix unitarization. For instance, 
the inelasticity reaches the asymptotic value 



2uo - 1 



as s 



(19) 



which is close to for oj close to 1/2. 

Our calculations for rj in (s,b) in the cases of the U- 
matrix and of the extended eikonal are shown in Fig. 
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FIG. 3: The allowed region (in white) for amplitudes x to be 
unitarised by an extended eikonal with ui = 0.525, together 
with the line corresponding to a hard pomeron amplitude with 
intercept 1.45. 




FIG. 4: Antishadowing effects in extended eikonal and U 
matrix schemes, for = 1.8 TeV (upper figure) and y/s — 
14 TeV (lower figure). 



representation for ui — 1/2. If, on the contrary, we let 
7 — > oo, we recover the extended form of the eikonal 
representation. Hence, when 7 varies from 1 to 00, and 
lu « 1/2, we obtain different forms of unitarization be- 
longing to a wide class with the similar asymptotic prop- 
erties, which we shall examine more closely in the next 
section. 



VI. COMPARISON OF THE BORN TERMS 

Another way to make the unitarization schemes lead 
to similar results is to use different inputs for z in (jHJ) or 
d2|). We can solve the equation 



-(1 - exp(iwx e ) = — 

U! 1 



Xu 



(21) 



to find which Born term in the U-matrix representation 
would give results similar to those of the eikonal. Writing 



Xu = Xu + iXu> and x, 



-x? 



+ ixi, we obtain 



2wxf = log 



tan(wxf) 



|x„|V 2 + 2xV + l 



| Xtt |2(u;-w') 2 -2 X £(w-w') + l 



(22) 



(23) 



This simplifies to a particularly simple expression in 
the case of a purely imaginary Xu- 



^xl = log 



xW 



1 



l- x i( w -w') 



(24) 



and the real part goes from to u>xf = 7r if (o> — tu')x u 
crosses 1. This relation is clearly discontinuous if u> 7^ u'. 
We illustrate this in Fig. 5 in the case ui = 1 (eikonal) 
and uj 1 = 1/2 (U matrix). At low energy, the phases 
are approximately the same. But at high energy when 
X^(s,6) — > 2, Xe has a discontinuity: its imaginary part 
goes to infinity, and its real part jumps by itt. On the 
other hand, the extended U matrix does not lead to such 
a singularity if u>' = 1. 



4. We see that both solutions have the same behavior in 
s and 6, but also that the extended eikonal has sharper 
anti-shadowing properties. 

Hence, the peculiar asymptotic properties of (|17[) are 
not unique to the [/-matrix pop . The extended eikonal 
(fTTj) has a similar asymptotic behavior for u> between 1/2 
and 1. 

Again, we can find a scheme that interpolates between 
these two forms Hal: 



G{s,b) 



[1 -(I- 



(20) 



If 7 



{l-iujx{s,b)h)~< 
1 this form coincides with the standard [/-matrix 



VII. NON-LINEAR EQUATIONS 

All previous schemes can be recast as non-linear equa- 
tions, which may be reminiscent of those obtained in 
QCD from gluon saturation. 

The simplest way [l9[ to get these is first to take the 
derivative of G with respect to x i n and (fT2")l : 



dG _ 1 
dx 



)G 



in the eikonal case, and 

dG 
dx 



(l + iuj'G) 2 



(25) 



(26) 
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xl 



FIG. 5: The relation between the Born terms of the U-matrix 
scheme and of the eikonal scheme, in the case of a purely 
imaginary \u- The real part of u\e is discontinuous and goes 
from to i-K at \ u = 2. 

for the U matrix. To make a tentative connection with 
saturation, we shall consider a purely imaginary Born 
term, and we shall write G = ig. Assuming \ = igox~ A , 
we can write the above equations as an evolution in y = 
log(l/x) at fixed b: 

^ = -log(l- Wfl )(l-w 5 ). (27) 
ay u> 

in the eikonal case and 

d f = A 3 (l - Jg) (28) 
dy 

for the U matrix. 

We see that the U matrix schemes lead to equations 
which look more natural than the corresponding ones in 
the eikonal case, as it is hard to imagine how saturation 
would lead to a log containing the amplitude. 

One can further generalize these equation to reproduce 
Eqs. (|15[) and (|20p . The corresponding nonlinear equa- 
tion will be 

%. = l±(l-(l- g yh)(l- U g). (29) 

ay uj 

It is easily seen that when 7 = 1 we obtain the non-linear 
equation for the U matrix. In the case u> = 1 Eq.(28) 
amounts to the standard logistic equation, and leads to 
the extended [/-matrix unitarization scheme. 



VIII. CONCLUSION 

In this paper we presented two new unitarization 
schemes which generalize the usual eikonal and U ma- 
trix unitarization schemes. We showed that they belong 
to two wide classes which cannot be mapped analytically 
one onto the other. We showed however that it is pos- 
sible to build a more general scheme which interpolates 
between the two. 

The basic behavior of G(s, b) as a function of s is 
mostly constrained by the value of to (or u>'), but not 



by the details of the unitarizing map. To illustrate this 
point, we show in Fig. 6 the behavior of G(s, b = 0) in our 
interpolating scheme (|20|) . for 7 = 9 (close to an eikonal) 
and for 7 = 1 (U matrix), for u> = 1 or 1/2. Clearly, the 
large- s behavior of the amplitude is controlled by uj, and 
not by 7. 




^ 5000 10000 15000 

O 2.h 




is (GeV) 

FIG. 6: The behavior of G(s,b = 0) from (gDJ for various 
choices of 7 and u. 

So the question of the asymptotic behavior of the elas- 
tic amplitude remains open. It is possible to build an 
infinite number of schemes in which the amplitude will 
saturate at the black disk limit, but there also exists an 
infinite number of schemes in which it will exceed it and 
eventually converge to the full unitarity limit. 

Up to now we do not have a decisive argument to 
choose one class of unitarization over the other. One pos- 
sibility would be to fit the existing data to determine 7 
and uj. It is however known that this is possible in eikonal 
schemes, so that it is unlikely that the constraints will be 
very stringent. However, the prediction for the total cross 
sections in these two classes of unitarization have large 
differences (see, for example [H|,[2(|) for the LHC energy 
region, and hence we may know soon which is realized. 
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